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REFLECTION FORMULA 
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ABSTRACT. The harmonic sawtooth map w (x) of the unit interval onto itself is defined. It is 
shown that its fixed points {x : w(x) = x} are enumerated by the n-th derivatives of a Meijer- 
G function and Lerch transcendent, serving as exponential and ordinary generating functions 
respectively, and involving the golden ratio in their parameters. The appropriately scaled Mellin 
transform of w (x) is an analytic continuation of the Riemann zeta function ¢(s) valid V—Re(s) ¢ 
N. The series expansion of the inverse scaling function associated to the Mellin transform of w (x) 
has coefficients enumerating the Large Schröder Numbers Sn, defined as the number of perfect 
matchings in a triangular grid of n squares and expressible as a hypergeometric function. A finite- 
sum approximation to ¢ (s) denoted by Cw (N; s) is examined and an associated function x (N; s) is 
found which solves the reflection formula Çw (N;1—s) = x (N; s) Ġw (N; s). The function x (N; s) 
is singular at s = 0 and the residue at this point changes sign from negative to positive between 
the values of N = 176 and N = 177. Some rather elegant graphs of the reflection functions 
x (N;s) are also provided. The Mellin and Laplace transforms of the individual component 
functions of the infinite sums and their roots are compared. The Gauss map h (x) is recalled so 
that its fixed points and Mellin transform can be contrasted to those of w(x). The geometric 


AE rel 
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counting function Nz, (x) = | of the fractal string Lw associated to the lengths 


of the harmonic sawtooth map components {wn (x)}°°_, happens to coincide with the counting 
function for the number of Pythagorean triangles of the form {(a,b,b + 1) : (b+ 1) < x}. The 
volume of the inner tubular neighborhood of the boundary of the map Lw with radius £ is shown 


dev(e)? —dev(e) 41 et oe 
2v(e) 2e à 
Also, the Minkowski content of Lw is shown to be Me, = 2 and the Minkowski dimension to be 
De, = 4 and thus not invertible. The geometric zeta function, which is the Mellin transform 
of the geometric counting function Nz, (x), is calculated and shown to have a rather unusual 
closed form involving a finite sum of Riemann zeta functions and binomioal coefficients. Some 


definitions from the theory of fractal strings and membranes are also recalled. 


to have the particuarly simple closed-form Vz, (€) = where v (£) = | 
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1. UNIT INTERVAL MAPPINGS 


1.1. The n-th Harmonic Sawtooth Function w,(z). 


1.1.1. Infinite Sum Decomposition. Let the n-th componentw,,(x) € [0,1]Vx € [0, 1] of the harmonic 
sawtooth function w(x) [9, 2.3] be defined as 


(1) Wns) =n(an+a—1)x (x, 1) 
where 
1 l <t ! 
ee ae 
0 otherwise 


is the characteristic function of the n-th harmonic interval(132). By setting n = |a~1| as in (134) 
we get the unit interval mapping 


w(x) =Wz- (£) 
= [z7] (æ [27t] +2- 1x (2,78) 
= am Wn(2) 
(3) =p nen + x — 1)x ( a) 


As can be seen in Figure 1, w(x) is discontinuous at a countably infinite set of points of Lebesgue 
measure zero 


H = fy: limy- w(x) 4 lim,,,+ w(x) 
(4) = fe ned} í ) 


The left and right limits at the discontinuous points are 


(5) 


limen- w(x) =1 
lim, e+ w(x) =0 


1.1.2. The Fixed Points of w(x) as an Iterated Function System. The iterates of the map 
(6) [w(@), w(w(x)), w(w(w(2))), w(w(ww(2)))), «J = [w (2), w%(2), w(x), we), z: 
have the form 


(7) w(x) =a, -— (x — c)br 
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FIGURE 1. The Harmonic Sawtooth map 


where {a,,b, € Z : r € N} is a pair of integer sequences and c € R is some constant. The sequence 
of quotients ¢* converges rather quickly to the fixed value 


(8) limpo F =x—c Va € [0,1] 
The explicit equation(sometimes called Schréder’s equation) for the fixed points of w(x), is 
Fix}, = {@: wnx) = 2} 


= {x : n(zn +z — 1)x(x, IŻ) =x} 
= {x :n(zn+xr-—1)= zr} 
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where G is the Meijer-G function and Fix,,(x) is the generation function 
Bie). So yo ikea 
(10) = ar 1 aa 


where ®(z,a,v) is the Lerch Transcendent (159), and ¢ is the Golden Ratio, which is the ratio 
of two numbers having the property that the ratio of the sum to the larger equals the ratio of 
the larger to the smaller. [7, Ch.XX][35, I.7][16, p.50][28] The number ¢ can be called the “most 
irrational number” because its continued fraction expansion, given by iterations of the Gauss map 
(78), converges more slowly than any other number. The constant ¢ satisfies the simple identities 


Bette a 
z ok l 
(12) g—$-1=0 


An interesting fact is that the density of a motif in a certain noncommutative space described in 
(23, 5.1] must necessarily be an element of the group Z + $Z. 


1.2. Integrals Transforms of w(x). 


1.2.1. Dirichlet Polynomial Series and the Mellin Transform of w(a). The Mellin transform of the 
harmonic saw map w(x), multiplied by 


stl 
s—1 
is an analytic continuation of the Riemann zeta function ¢(s)V — Re(s) ¢ N. This form of the 


zeta function, denoted by ¢.,(s), is the infinite sum of the Mellin transformations of the component 
functions. 





(13) t(s)=s 


A Wn(x)a*—lda 
n+l wp (a7 a7" "der 
= Jonen +z -— 1)x(z, IŻ”) tdg 
(14) =f", _n(en +z- 1)r® tdr 
daoa- ne elat) ael) 


s? +s 
— nén+s(n+1)>—(n+1)in 











= ns(n+1)*-+n%s*(n+1)° 
n°? —n(n+1)7%—sn7* 
s2ts 


There is a conjugate pair of inverse branches of 7(s 


) 
Te (t) ={s:7(s) =t} 


found by solving 








(15) = s: 9S} =F 
t 1 viote 
27 DS 2 





where 7} '(t) and TZ} (t) denote the positive and negative solutions respectively. The coefficients 
in the series expansions are integers enumerating the large Schröder numbers Sn which count the 
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number of perfect matchings in a triangular grid of n squares, named after Ernst Schréder(1841- 
1902).[17, A006318][4, p.340][10] 


(16) Sn =n (37! Pom aa ja 


where pF is a hypergeometric function. We have 





nt) 0 n=0 
limo H = -1 n=1 
—Sn n>2 
(17) nay, =I n=0 
lim;_59 Calis ®© = 2 n = 
Sn n22 
(18) 
n |0|112]13] 41516] 7 8 9 10 11 12 13 14 | 





















































Sn | 1} 1} 2} 6] 22 | 90 | 394 | 1806 | 8558 | 41586 | 206098 | 1037718 | 5293446 | 27297738 142078746 | 





The residue at the singular point s = 1 of r(s) is 





Res (7(s)M[wpn(x); £ > s]) = tes 


s=1 s 


s? +s n! —n(n+1) sn? 
—1 s?@+s 


n(n+1) -nt +sn R ) 





(19) 























Cw (s) = ¢(s) 
= 1(s)M[w(x);2 > s] 
=1(s) f) w(x)a>-!de 
= 9281 Ji [a3] (a [2-1] +e- 1) at ae 
(20) = 8th Ppa M[wr(z); £ — 5] 
1 
= sft Eri J a nan +a ije de 
nPI 
oo P E eee 
œo  n(n4+1) =n! psn7® 
E 2 n=1 ! ? s—1 
= = yeaa n(n + 1)7§ —n'-§ + sn75 





Unlike the Mellin transform of the Gauss map (78), which must be multiplied by the factor s then 
subtracted from -4 before it equals ¢(s), the “harmonic sawtooth continuation” ¢,,(s) of the zeta 
function ¢(s) has the fortuitous property that it equals ¢(s) after multiplying M|w(x);2 —> s] by 
T(s). This property of w(x) allows us to put 7(s) inside the sum to get an expression, denoted 
by ¢w(N; s), involving the difference of two Dirichlet polynomials, one of which is scaled by s. The 
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substition oo — N is made in the infinite sum appearing the expression for ¢,,(s) to get 


Cw(N;s) =7(s) pean M [wn (x); x > s] 
= Ay Pi hA n + sn“ 
œa) = Ay (s+ E 5D nt — EN ns) 


= N cos(7s)V(s—1,N+1) , * 
= GINI T(s) + C(s)V¥s EN 








with equality in the limit except at the negative integers 


(22) limy soo Gw(N3s) = ¢(s)V — s g N* 

The functions Çw (N; s) have real roots at s = 0 and s = —1. That is 

(23) lim ¢, (N; s) = lim ¢, (N;s) =0 
s—0 sol 


The residue of ¢,,(V;s) at s = 1 is given by 


Res(Gw(Nis)) = Res (7(s) Dns Mlwn(e):2 > 5) 


(24) = Dru Res (roluna: z> 4) 


= Sia ais 
= n=1 n2+n 
N 





N+1 
Thus, as required 
limno Res(Gw(Nis)) = limysoc Ye we 

25 ; 
(25) = limy oo x5 

=1 
The function 7(s) has zeros at —1 and 0 and a simple pole at s = 1 with residue 

= s+1\ _ 

(26) Res (r(s)) = Res (st#t) =2 


The Mellin transform of r(s) has an interesting Laurent series, convergent on the unit disc, given 
by 


M[r(s);s >t] = t(s)s'~'ds 


= [°° gstigt-1 
= fo 8 ss ds 


= a 2n—3 

en = Salam ye Vel 
= 4i nE aie: ee ) 

cos(7t) 


sin(7t) 








= —2r 


The transformations M[w,,(a);2 — s| have removable singularities at —1 and 0 where the limits 
are given by 





Jim M[wn(2);a > s] =n? In ) + nin (24+) -n 


(28) 1 " 


Ger 
+ nln (#4 


lim M [wn (2); — s] = 
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lity 35 e~ Ml[wn(z);x—s] = (n+ 1)?n "e7! 
2 


(29) lim, +_1 e Mlwn (x);2— 8] = n” +n) (n 4 1)? en 


So, (14) can be rewritten as 

















n? In (2+) +nIn(2*)-n s=-1 
(30) M [wns]; s>s] ={ 1tnm(SP) s=0 
z zenti — otherwise 
Furthermore, we have the limits 
limno lim M [wn(£); £ > s| = limpo lims+-1 aen oa 
so 
(31) = limps n? In (==) + nin (==) n 
=} 
and 
i : _ or ' n! —n(n+1) -sn ° 
limy—soo lim M [wn(x£); £ > s| = limno lims+4o aa 
s> 
(32) = limps 1 + nln (24) 
=0 
The sum of limits over n at s = —1 is Euler’s constant. [13][14, 1.1] 
co)=60M [wn (x);235-1 co l+nln z 
nal Mania a n= intel 
= lim, +1¢(s) — 4 
(33) = limn Dee £ — In(n) 


=) 
= 0.577215664901533... 


1.2.2. The Reflection Formula for Gw(N;s). There is a reflection equation for the finite-sum ap- 


proximation ¢,,(N;s) which is similiar to the well-known formula ¢(1—s) = x(s)¢(s) with 
x (s) = 2 (2r) * cos (Z£) T (s). The solution to 
(34) Cw (N; 1 = 8) = x (N; 8) Cw (N; 8) 


is given by the expression 


w(N;1— 
x (N58) ANT 
S^, -n+(n41)° nt 


nsl ns- ls 








(35) = N Sal epa) nin ss 





n=1 3-1 
(s—1) N_, -n° +(n+1)° tnnt —n: s 
s > N_, ants 4(n+1) Snn ss 








which satisfies 

(0) x(N;1 = 8s) =x(N;s) t 

The functions y (N; s), indexed by N, have singularities at s = 0. Let 
a(N) = Di n(n (n+ 1) — ln (n)) 

(37) b(N) = D a e 


c(N) =1~% n (In (n+1)?—In (n)”) 
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then the residue at the singular point s = 0 is given by the expression 
(38) 
; = .g)-1 
Res(x(N;s)) = —Res(x(N;s)~) 
N(n(P(N+1))=cUW)) 
1t+7+W(n+2)— wr +0(N) (WV =a) (NF) 
> a(N)- N 
1 N 2 2 
Ly © (n+2)— p2 +N, Pn? antn? nln) N (WTNH) -3 Dh n(n? ~In(m)?)) 
N+1 n=1 n(n+1) (W=SA_, n(n(n+1)—In(n)))(W+1) 

n ( Nai n(In(n-+1)—In(n)))—N 























which has the limit 


(39) lim Res(x(N;s)) =1 


N-oo s=0 
We also have the residue of the reciprocal at s = 2 


2N =- —2W(1,N+1)+2¢(2) 
40 Res(x(N;s)-") = way" 
( ) s=2 (x ) ) S4 NLE n(In(n41)41n(n+1)n—In(n)—n In(n)) 








which vanishes as N tends to infinity 


(41) lim Res(x(N;s)~') = 0 


N= s=2 


As can be seen in the figures below, the residue at s = 0 changes sign from negative to positive 
between the values of N = 176 and N = 177. 











FIGURE 2. {Res(x(N; s)):N= 1...250} 
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FIGURE 3. {Res(x(N; s))-1:N=1... 250} 


For any positive integer N, we have the limits 
lims_+0 x (N; s) =0o 
lims0 gX (N38) = 00 
ua) impala) = 
lims—1 x (N; 8) =0 
lim,+2 x (N; s8) =0 
lims41 &x(N;s) =0 


The line Re (s) = 4 has a constant modulus 























1 
(43) x (wig tis) [= 
There is also the complex conjugate symmetry 
(44) x (N; x + ty) = x (N; z — ty) 
If s =n € N* is a positive integer then x (N; n) can be written as 
: — w(N;l-n) 
x(N;n) = teo NG) 
= er ee a) 
(45) = e SAE ALNET +C(n) 
_ Ena (n-m +m” -(m+1)” m) 
a cosent mN ED 4 ¢(n) 


The Bernoulli numbers[1] make an appearance since 


(46) X(N; 2n) Gu(N32n) = Ban (N +1)? ZH +... 


250 
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The denominator of x (N; n) has the limits 


(47) limy Cw Vin) =¢(n) 


mnoo Gu (Nin) =1 


Another interesting formula gives the limit at s = 1 of the quotient of successive functions 























lim x(N+1;s) = (N+2)N(N+1—a(N+1)) 
s=1  X(N;s) ee a(N)) 
(48) O (N+2)N(N41 eS) 
2 


(N+1)?(N- n(In(n-+1)—In(n))) 






































































































































Let 
(49) v (s) = x (00; 8) = 


Then the residue at the even negative integers is 


(50) Res (v(s)) = 


ca-n) 
zaag- reven 
1.2.3. The Laplace Transforms L|wn(x); x —> s]. The Laplace transform L[wn(£);x —> s] and its 
roots are calculated to shed light on the behaviour roots of the Mellin transforms M|wn (x); £ — s] 
but it is unclear whether this is accomplished. The Laplace transform (135) of the n-th component 
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FIGURE 5. {x(N;s):s=4...2,N = 1... 100} 


function is given by 


(51) 


Liw); £ >s] = = Jo ae 73 dy 
= = fén rn ' x — 1)x(x, IŻ )e-**dx 
= =f n(an + x — l)e "sdg 


_ niate” mT —(n?4+n+s)e~ a 
= 2 





There is a removable singularity at s = 0 which has the limit 


L[w,(x);2 > 0] = lims_o L[w e Ea z> s] 














n(n+1)e a 
(52) = lim ae 
= zRes (T Mwn sa > s]) 
= mnt 
Additionally, 
al L[wn (a); cs 0] Z Da lim,_s0 n(n es tn+s)e n 
Co 
= Anai ea 
(53) = limn—oo m Mlwn (x); x > s] 


ni~$—n(n+1)7%—sn75 
s?+s 





= limy +o lims 4-1 
1 


2 
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The roots of L[w,(x);x — s| are enumerated by 
(54) pk (m) ={s: L[w,(z);2 > s| = 0} = -n(n + 1)W(m, —e7") — n(n + 1) 
where W (m, x) is the Lambert W function (146) and m € Z,n € N. It can be verified that 


emt DFW (m-e DW (m,—e71) perm CHW (om, =e!) 


n(1+W(m,—e-"))2(n41) =0 





(55) Ljwn(z);£ => p% (m)| = 
where —e~! is expressed as the continued fraction via its quotient sequence 
(56) Seo) = |-1,1,1,1,2,1,1,4,1,1,6,1,1,8,1,1,10,1,1,12,1,1,14,1,1,16,1,1,18,1,1,20,...] 


The roots p% (m) satisfy a functional reflection equation with respect to m 





-n(n + 1)W(m,—e7!)— n(n +1) =-n(n+1)W(—-m — 1,-e7!) — n(n + 1) 
where = Re(x) — Im(x) denotes complex conjugation. The quotients of the roots of consecutive 
transforms is 








Pron (m) n 
(58) an) T TI 
Thus 
: Prion (M) . 1 
(59) limn—oo BE, (mn) = limnsoo Fj = 1 


1.2.4. The Roots pM (m) of M[w,(x);x > s]. Define 


M,y(8,7) == T(s)M[wn(x); x => | 
(60) = n(n+1)~8§—n'~S$+sn7* 





as in (21) and its infinte number of inverse branches (which are currently lacking closed-form 

expression if such a thing is possible except when n = 1), where the branches are indexed by m 
Mz (z,n,m) ={s:7(s)M[wn(x);2 > s| = z} 

n(n+1) =n! sn? 


s—1 


(61) 





= g = z 
then we see that the first function where n = 1, M7} (z,1,m), has the closed-form 


W (m, BOP ) +1n(2) 














(62) Myu (z, 1, m) = In(2) 
where W is the Lambert W function (146). It is verified that 
Ww (m BYD) tme) w (m, BOD) 42 
Ži E ae TH) = 
(63) Mrw(Ma5(z,1,m),1,m) = TEA =ž 





mO) 
Furthermore, let p% (m) denote the m-th root of M[wn(x); £ > s] 


: M[wn(£); x > s] = 0} 


(64) Pir me 7 oe n, m) 


which satisfies 
(65) Im (pi! (m)) > Im Ce (m)) 
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iP, _[m]:n=1..9,m=-00..-1} 
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0 50 — " eo Es = 200 250 
FIGURE 6. The Roots p% (m) of {L[wn(z);2 > s] :n€1...9} 
(66) limm+too Re(p! (m)) =0 
Thus, 
(67) limmo arg (ont. (m)) =f 
1.2.5. Quotients and Differences of pi, (m). Let 
(68) Apu, (M) = py, (m+ 1) = pur, (M) 


be the forward difference of consecutive roots of M|wn(£); x — s]. The limiting difference between 
consecutive roots is the countably infinite set of solutions to the equation n? + (n +1)? = 0 given 
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by 








+ 
p (m) 
= Mmo €m 


2 277% 
~~ M[x(z,I¥);x—0] 
27% 





ween) 


s 


lims>o( 
2ri 
In(n+1)—In(n) 


Let Qa denote the limit 


wn 











pat, 5 limmo pa Gm) 
Apis, (+00) 
(70) = Ar, 1 (oo) 





mix(e(2 iq) 0) 


In(n)—In(n— 1) 
In(n+1)—In(n) 


then we also have the limit of the limits Qis as n — oo given by 


wn 














In(n+1)—In(n) 


i g Ap™ (+00) 
limn—++oo Qu — limno aM a 
p Ap (+00) 
wn wn E 
In(n)—In(n—1) 
as ir HEDDET 
(71) = n—+too | mmn-D-mnmn-3 
“In(ay—In(n—3) 
Į -l —1 
= Tir, duh (SSS) 








= 
—pM (m 
The limiting quotients £ a ne = Pin (™ Pin MI) ag m — o are given by 
(72) 
í M (m)—-oM (m 7 = a g p ; p 
liMmm—oo ePwnl )=Pwn (MEI) =1-—2sin (eonim) 2i cos (tm) sin (son) 


Qri 


=e mmFi)-In(n) 


where we have 


: M „M ; M _,M 
| liMm—oo eP¥n (m)=Pwn Pn) = lima |ePun (™)— Pain (m+1)| 


(73) = Ve tee 


=1 
and 


(74) limpo litt oq elun -Pun MH) — 1 


600 —f 
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FIGURE 7. {pM (m):m=1...5} 
1)M[wn(ax);x — s] defined by 


1.2.6. The Laplace Transform L[|(s— 1)M[wn(x£);x > s];s > t]. The Laplace transform of (s — 
Lj(s —1)M[w, (x); x > s]; s — t] 
(75) 


= Lin(n +1) + sn™° — nts; s — t] 
=f 
t 


1 
o n(n +1) pen — ntse ds 
has poles at — ln(n) and — ln(n + 1) with residues 








Res 
(76) 


| In( oie )t+In(n +1)+n In(n)?—n In(n) In(n+1) 
(n(n 
t=— In(n) 


(n)+t)? (In(n+1)+) 





(Lis — 1)M [wn (x); £ > s]; s > tl) 


er (Lis —1)M|[w,(2);2 > s]; s > tl) =n 


1.3. The Gauss Map h(x). 


(77) 


formation, maps unit intervals onto unit intervals and by iteration gives the continued fraction 


1.3.1. Continued Fractions. The Gauss map h(x), also known as the Gauss function or Gauss trans- 
expansion of a real number [38, A.1.7][35, I.1][12, X] The n-th component function h,,(x) of the 
map h(x) is given by 


hn(2) 


— i-xn 





x (x, I3) 
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The infinite sum of the component functions is the Gauss map 


A(x) = p= hn(2) 
(78) = Dna "x (2, Ef) 


I 


02 03 04 05 06 07 0.8 0 


= hO 





FIGURE 8. The Gauss Map 


The fixed points of h(x) are the (positive) solutions to the equation hn(x) = x given by 
Fix, = {x : h(x) =a} 
(79) 
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1.3.2. The Mellin Transform of h(x). The Melin transform (138) of the Gauss map h(x) over the 
unit interval, scaled by s then subtracted from $, is an analytic continuation of ¢(s), denoted by 
¢n(s), valid for all— Re(s)) g N. The transfer oaii and thermodynamic aspects of the Gauss 
map are discussed in [42][41][40][39][36]. The Mellin transform of the n-th component function 
Wn (x) is given by 


M[hn(x);¢ > 8] = fy An(x)a*-lde 
= } tomy (x, I7) x5 ldg 
(80) = ee (£7! — [æ71]) ada 





ll 
=< 
er 


l-—an,,s—1 
=r dr 


n(n+1)-°+s(n+1)-—n!-: 


s“—s 





which provides an analytic continuation Ca (s) = ¢(s)V(— Re(s)) ¢ N 
a(s) = 





se? a] 
-sfo h(x) i 


= sy -$ fg (01 = [2 1) ode 
(81) = aoit up (z); z > s] 


a at =T n= 1 (n(n ie F a(n A? ni) 


i l 














= Gen Ge si n=1” 
1.4. The Harmonic Sawtooth Map w(x) as an Ordinary Fractal String. 
1.4.1. Definition and Length. Let 
(82) H = (4,3) 


be the n-th harmonic interval, then {w(x) € Ly : x € Q} is the piecewise monotone mapping of 
the unit interval onto itself. The fractal string Lw associated with w(x) is the set of connected 
component functions w,(x) C w(x) where each w,(x) maps// onto (0,1) and vanishes when 
x ¢I#. Thus, the disjoint union of the connected components of Lw is the infinite sum w(x) = 
S>] Wn(x) where only 1 of the wn (x) is nonzero for each x, thus w(x) maps entire unit interval onto 
itself unigely except except for the points of discontinuity on the boundary OL, = {0, + :n€N*} 
where a choice is to be made between 0 and 1 depending on the direction in which the limit is 
approached. Let 


(83) W(x) =n(an+a—1)x(a, IŻ) 


where x(x, IŻ) is the n-th harmonic interval indicator (127) 





(84) xæ) = 0 (st) — 0 (a=) 


The substitution n + |ż| can be made in (132) where it is seen that 


(85) x(x tf.) =o (H) a(i) =: 
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and so making the same substitution in (83) gives 


l 


02 03 04 05 06 07 08 0.9 


= ea 


(86) 


FIGURE 9. The Harmonic Sawtooth Map 


The intervals I” will be defined such that wy, = |I?| = |wn (x)|. Let 


(87) E ES 





INTEGRAL TRANSFORMS 19 


be the midpoint of IŽ then 











88 Te = (tn — gm bg + Mayet) 
(88) 4n+1 4n+3 
4n(n+1)? 4n(n+1) 
so that 
lwn = |Wn(x)| 
1 
= joe +a2—1)x(x,I#)dzx 
a ES w(x)dz 
(89) =|] 
1 w 
= fo X(x, Ip de 
— _4n A4n+1 
4n(n+1) 4n(n+1) 
~~ 2n(n+1) 
20000 
15000 
1 
Fw 10000 
5000 





: 7 : | 
10 20 30 40 50 60 70 80 90 100 
n=1..100 


FIGURE 10. Reciprocal lengths ¢w; + 


The total length of Lw is 


Lu] = fo w(z)dz 
(90) = age Mh 
= X n=1 re) 
~2 
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FIGURE 11. x(a, I") green and w,(x) blue for n = 1...3 and x = 5.1 


1.4.2. Geometry and Volume of the Inner Tubular Neighborhood. The geometric counting function 
(109) of Lw is 


Nc (2) = #{n 2 1: lwz! < a} 


w 


(91) - [aa T +1)n <a} 


which is used to calculate the limiting constant (111) Cù appearing in the equation for the Minkows- 
ki content 








(92) =I 


The function Nz, (x) happens to coincide with [17, A095861], which is the number of primitive 
Pythagorean triangles of the form {(a,b,b +1): (b+1) < n}. [6, 171-176][37, 10.1][18, 11.2-11.5] 
Let 


(93) v(e) = min(j : fw; <2e) = || 
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which is the floor of the solution to the inverse length equation 


Etv ette YZE = {n : lwn- = 2e} = {n : yhy = 2e) 
á etve? +e HF 


E — 
2 =€ 


(94) 


Then the volume of the inner tubular neighborhood of 0£,, with radius € (108) is 














Vele) = 2€Ne. (4) + Dy bw 
Co 1 
= 2eNz,, (+) T > n=v(e) 2(n+1)n 
a = 2eNc, (ze) + we 
= 2e at 1 EO 
— 4ev(e)? —4ev(e)+1 
5 2v(e 
since 
(96) ea CESA = Ia 
and by defintion we have 
lim,_,9+ Ve, (€) =0 
(97) lime 00 Vee (€) = |Le] = } 


Thus, using (92) and (95), the Minkowski content (110) of Lw is 

















: V£w (E 
Mge, = lime;o AO] 
l 1 2 Viti 1], 1 | etWerte | 
= HMe_,o+ Te | 4€ 2 2| 7 2 2e 
(98) — Cy? tw 
va) 1—Dew 
21-5 
peas 2 
127 
2 


1.4.3. The Geometric Zeta Function Cc (s). The geometric zeta function (112) of w(x) is the 
Dirichlet series of the lengths @w,, (89) and also an integral over the geometric length counting 
function (91) Ng, (x) 


eu) SEL 
= Ss (zem) 
(99) =) p12 (etl) n: 
=s fj Ne, (x)a7 tde 
2 fo° [AEE — a 
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FIGURE 12. Geometric Counting Function Nz, (x) of w(x) 


The residue (115) of ¢¢,,(s) at De, is 








Res(¢r(s)) =lim, pš (s — Deels) 
s=De ad 
= lim, „+ (s — 4) ce(s) 
(100) = lim, 1+ (8— 3) Dnc 2 (n + 1) en 
=lim, 34 ben |e = 4] ode 


The values of ¢¢,, (n) at positive integer values n € N* are given explicitly by a rather unwieldy sum 
of binomial coefficients and the Riemann zeta function ¢(n) at even integer values. First, define 


(n—1)(1—(-1)"* 











an = 5 
—1)” t1 (n—1 1)” 
(101) Dn — (=) ms Die 4 CD 
cn = (—-1)"(n-1) 
dn = = 


then 


2 
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min(n: tw <2e) 





(0) 0.05 0.10 0.15 0.20 0.25 
TETE S 
1000” 4 
FIGURE 13. {v(e) = min(n : lwp < 22) : € = io- Gh 
—1)7 (2271 2(—1)" ("ten -dn+3)e(dn 42n- 3—-2M-—cn 
(102) Ce (n) = (1) Ca) + pDr Cyr Emaar 5 n=8-2m—cn) 





a 
The terms of c (n) from n = 1 to 10 are shown below in Table 1. 




















1 
2 
-= Re 
to aO À 
Ig +3 ¢ (2) +g C (4) 
+R -ECD -peA 
-Al + $3 ¢ (2) +35 ¢ (4) +56 (6) 
+% Brea) -Bou)  - £608) 
-sbs iea +H) 4 866) +c) 
-OE e E po 
Bio Aco ibea elo +E elet d cao) 
TABLE 1. {¢c,,(n) = ÈE rown:n=1..10} 


24 CROWLEY 


0.5 
0.4 
0.3 
VL © 
w 
0.2 
0.1 
0 0.02 0.04 0.06 0.08 0.10 0.12 
=0..— 
aa 


FIGURE 14. Volume of the inner tubular neighborhood of 0£,, with radius e {Vz,,(€):¢=0...§} 


2. FRACTAL STRINGS AND DYNAMICAL ZETA FUNCTIONS 


2.1. Fractal Strings. A a fractal string £ is defined as a nonempty bounded open subset of the 
real line £ C R consisting of a countable disjoint union of open intervals I; 


(103) £= (5 
j=l 
The length of the j-th interval J; is denoted by 
(104) L; = |G 
where |-| is the 1-dimensional Lebesgue measure. The lengths @; must form a nonnegative monot- 


ically decreasing sequence and the total length must be finite, that is 


[Lh = Ej b < 00 
(108) hb... 2 Ghj O0 


The case when £; = 0 for any j will be excluded here since £j is a finite sequence. The fractal string 
is defined completely by its sequence of lengths so it can be denoted 


(106) fate, 
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2 
1.9 
1.8 
VL (e) 
— =| AF 
Je 
1.6 
1.5 
0 0.02 0.04 0.06 0.08 0.10 0.12 
e=0 = 
-3 


FIGURE 15. [l :e=0...3} and Ee = 2 


The boundary of £ in R, denoted by ðL C Q, is a totally disconnected bounded perfect subset 
which can be represented as a string of finite length, and generally any compact subset of R, also 
has this property. The boundary ðL is said to be perfect since it is closed and each of its points 
is a limit point. Since the Cantor-Bendixon lemma states that there exists a perfect set P C OL 
such that ðL — P is a most countable, we can define £ as the complenent of OL in its closed convex 
hull. The connected components of the bounded open set £\0L are the intervals J;. [25, 1.2][32, 
2.2 Ex17][23, 3.1][30][21][20] [15][22][11][19] [29] 


2.1.1. The Minkowski Dimension De and Content Mz. The Minkowski dimension Ds € [0,1], 
also known as the box dimension, is maximum value of V (e) 


(107) De =inf{az0:V(e)=O(e “*) ase => 0t} = ¢e(1)= X} £ 


n=l" 
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where V (e) is the volume of the inner tubular neighborhoods of ôL with radius € 


Ve(e) = |r € L£: d(x, OL) < el 


Lj 22e £3 <2e 
= y 2e + J j 
J J 
Lj <2e 


1 


and N(x) is the geometric counting function which is the number of components with their recip- 
rocal length being less than or equal to zx. 


(108) 


(109) 


The Minkowski content of £ is then defined as 











Me =lim. 49+ 3562 
Crp2i-Pc 
(110) = 
— Res(ce(s);De)2 Pf 
oa Dce(i-De) 
where Cz is the constant 
(111) Ce = Timg yoo “$2 


If Mc € (0,00) exists then £ is said to be Minkowski measurable which necessarily means that the 
geometry of £ does not oscillate and vice versa. [27, 1] [3][24][26, 6.2] 


2.1.2. The Geometric Zeta Function ¢c(s). The geometric Zeta function ¢c(s) of £ is the Dirichlet 
series 


Ce(s) = don G 
(112) =s T Nea dy 


which is holomorphic for Re(s) > Dg. If £ is Minkowski measurable then 0 < Dg < 1 is the simple 
unique pole of ¢¢(s) on the vertical line Re(s) = Dg. Assuming Çz (s) has a meromorphic extension 
to a neighboorhood of Dg then ¢c(s) has a simple pole at ¢e (De) if 


(113) Ng (s) = O(s?£) ass > 00 

or if the volume of the tubular neighborhoods satisfies 

(114) Vele) = O(et7P£) ase = OF 

It can be possible that the residue of ¢-(s) at s = Dg is positive and finite 
(115) 0< p a — De)éc(s) < 00 


even if Nc(s) is not of order s?£ as s > œo and Vz (e) is not of order ¢1~?£, however this does not 
contradict the Minkowski measurability of £. 
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2.1.3. Complex Dimensions, Screens and Windows. The set of visible complex dimensions of £, 
denoted by Dc(W), is a discrete subset of © consisting of the poles of {¢c(s) :s € W}. 


(116) De(W) ={weW: ¢c(w)isapole} 


When W is the entire complex plane then the set De (C) = Dg is simply called the set of complex 
dimensions of £L. The presence of oscillations in V(e) implies the presence of imaginary complex 
dimensions with Re(-) = Dg and vice versa. More generally, the complex dimensions of a fractal 
string £ describe its geometric and spectral oscillations. 


2.1.4. Frequencies of Fractal Strings and Spectral Zeta Functions. The eigenvalues Apn of the Dirich- 
2 

let Laplacian Au(x) = — u(x) on a bounded open set Q C R correspond to the normalized fre- 

quencies f, = ¥ an of a fractal string. The frequencies of the unit interval are the natural numbers 


n € N* and the frequencies of an interval of length £ are né—'. The frequencies of £ are the numbers 





(117) frj = he; WV b,j e N* 
The spectral counting function N,,¢(x) counts the frequencies of £ with multiplicity 
(118) Aaaa = j= Ne (3) 

T j= [z£] 


The spectral zeta function ¢,¢(s) of £ is connected to the Riemann zeta function (??) by 


Gue(s) = Dogar Deja hE 
(119) = €(8) nai G 
= ¢(s)¢c(s) 
(27, 1.1][25, 1.2.1] 


2.1.5. Generalized Fractal Strings and Dirichlet Integrals. A generalized fractal string is a positive 
or complex or local measure n(x) on (0,00) such that 


(120) ie n(a)da = 0 


for some zo > 0. A local positive measure is a standard positive Borel measure 7(J) on (0,00) 
where J is the set of all bounded subintervals of (0,00) in which case n(x) = |7(a)|. More generally, 
a meausre 7(x) is a local complex measure if n(A) is well-defined for any subset A C [a,b] where 
[a,b] C [0,00] is a bounded subset of the positive half-line (0,00) and the restriction of 7 to the 
Borel subsets of [a,b] is a complex measure on fa, b] in the traditional sense. The geometric counting 
function of 7(a) is defined as 


(121) Ny(x) = fê n(a)de 
The dimension D, is the abscissa of conergence of the Dirichlet integral 
(122) Gaile) = | 2-eIn(a)laz 


In other terms, it is the smallest real positive o such that the improper Riemann-Lebesgue converges 
to a finite value. The geometric zeta function is defined as the Mellin transform 


(123) G(s) = f anade 
where Re(s) > Dy. 
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2.2. Fractal Membranes and Spectral Partitions. 


2.2.1. Complex Dimensions of Dynamical Zeta Functions. The fractal membrane Tz associated 
with £ is the adelic product 


(124) =F 


where each 7; is an interval J; of length log(£;*)~*. To each 7; is associated a Hilbert space 
H; = L?(I;) of square integrable functions on Ij. The spectral partition function Zc(s) of £ is a 
Euler product expansion which has no zeros or poles in Re(s) > Dm(£). 
Zc(s) = ie TE 

= 15 ae Ze; (s) 
where Dm(£) is the Minkowski dimension of £ and Zg (s) = z4x is the j-th Euler factor, the 
partition function of the j-th component of the fractal membrane. (23, 3.2.2] 


(125) 





2.2.2. Dynamical Zeta Functions of Fractal Membranes. The dynamical zeta function of a fractal 
membrane £ is the negative of the logarithmic derivative of the Zeta function associated with £. 
Ze(s) = -8 In(Ce(s)) 


E ce(s) 


> Ce (s) 


(126) 


3. SPECIAL FUNCTIONS, DEFINITIONS, AND CONVENTIONS 
3.1. Special Functions. 


3.1.1. The Interval Indicator (Characteristic) Function x(x, I). The (left-open, right-closed) inter- 
val indicator function is y(x, I) where I = (a, b] 


1 «el 
(127) _fla<a<b 
-~ ) O otherwise 





= (x — a) — O(a — a)0(ax — b) 
and @ is the Heaviside unit step function, the derivative of which is the Dirac delta function 6 
[stelde- =O(2) 
(128) _f 0 «2<0 
A(z) ~ ) 1 «#30 
The discontinous point of 0(x) has the limiting values 
lim,_.9- O(x) =0 
co?) lim, 0+ A(t) =1 


thus the values of x(x, (a,b)) on the boundary can be chosen according to which side the limit is 
regarded as being approached from. 


Ee i = 6(a — b) 
(130) i Se (a ee 
lim, 5+ x(a, (a, b]) = 
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3.1.2. “Harmonic” Intervals. Let the n-th harmonic (left-open, right-closed) interval be defined as 
(131) Ñ = (t) 


then its characteristic function is 








6 
(132) =o 
1 

2 ESI A 

0 otherwise 
As can be seen 
Co = UR (A = (0,1) 
Texel) = Drax (e (zl) =x, Ou) 


The substitution n + |ż| can be made in (132) where it is seen that 


(133) 


I 





H p slet] +x—1 z|e~*|-1 = = 
(134) x (cyte) =0( pa ) o( -= ) =1 Yze[-1,+1] 
3.1.3. The Laplace Transform L?[f (x); x — s]. The Laplace transform [2, 1.5] is defined as 


(135) L?[f(a);a > s] = f? f(w)en*da 


where the unilateral Laplace transform is over the interval (a,b) = (0, 00) and the bilateral transform 
is over (a,b) = (—co,00). When (a,b) is not specified, it is assumed to range over the support of 
f(x) if the support is an interval. If the support of f(x) is not an interval then (a,b) must be 
specified. Applying L to the interval indicator function (127) gives 


L? [x(x, (a,b));£ > s] = f? x(x, (a,b))e~*dax 
f (O(a — a) — O(a — b)0(a — a))e~*dax 





II 








i a 
Pets etje 7+0) 
The limit at the singular point s = 0 is 
lims so LP [x(e, (a, 6));2 > s] = lim, 49 — ee” 
(137) a ; ) 
=O0—-4 


3.1.4. The Mellin Transform ML[f(x);x — s]. The Mellin transform [31, 3.2][33, II.10.8][5, 3.6] is 
defined as 


(138) Mf (a);2 > s] =f) f(w)ae de 


where the standard Mellin transform is over the interval (a,b) = (0,00). Again, as with the notation 
for Laplace transform, the integral is over the support of f(a) if the support is an interval and (a, b) 
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is not specified, otherwise (a,b) must be specified. Applying M to the interval indicator function 
(127) gives 


Mè[x(x, (a,b));2 > s] = L x(x, (a, b))as~tdax 
(139) = f (0(x —a) — O(a — b)O(x — a)) asda 








The limit at the singular point s = 0 is 


Mz (x(a, (a, 6));2 > 0] = M [x(x, (a, b));æ > 0] 
= limo MELy (x, (a, b)); £ > s] 


5—a? 


(140) 





= lims+o0 


= In(b) — In(a) 
The Mellin transform has several identities [31, 3.1.2], including but not limited to 


M|f (az); x > s] A RO. 
M|x* f(x); x£ > s] Mf REE 
M[f(x*); £ > s] = gif f(z); > §] 

(141) M[f (£72); x > s] = -M | f(z); Aaa 
M{a f(z”); x — s] = 7 f(z); xz > ie 
Miztf(e*);2 > 5] = 1M |f(x)2 > -e| 
M 





ln(x)” f(z); >s] = 2M [fa z); — s] 
where a > 0, u > 0, and n € N. The Mellin transform of the harmonic interval indicator function 
(132) is 








M [x (s,f); > 0] = S3 x (e, (z $) ) ae tae 
(142) = J7 (8 (t-z) - 9-2) ear 
a 


which has the limit 
M Ix (x, TF) £> 0] = lim, 9 M Ix (x, TF) 3m s] 
(143) = lim, 45 E0 
= In(n + 1) — ln(n) 
The Mellin and bilaterial Laplace transforms are related by the change of variables x > — In(y) 
resulting in the identity [31, 3.1.1] 
MPL- a jas] = LIS i ys 


(144) Ff f(—In(a \as—ldar Zs y)e ~¥S dy 


3.1.5. The Lambert W Function W (k,x). The Lambert W function [8][34] is the inverse of xe” 
given by 
W(z) ={a: xe” = z} 
= W(0, z) 
= 1+ (In(z) - 1 exp ($ fo” shy m (SERS Jaz) 
k-1l yk 


co =k z 
a oat ae 


(145) 








INTEGRAL TRANSFORMS 
where W (a, z)Va € Z, z € {0, —e7}} is 
(146)  W(a,z) =1+ (ina + In(z) — 1) exp (5 fo? yn (BGS eee) ae) 
A generaliztion of (145) is solved by 


(147) {r:r =z} = a 











The W function satisifes several identities 
W(zjeWO =z 


W(zIn(z)) = ln(z) Vz <1 
yar = A 
EMANE = ZW z)” 

(148) In(W(n, z)) = ln(z) — W (n, z) + 2inn 
w (-22) = — ln(z) Yz € [0, e] 
womg) yt 


—In(z) 


where n € Z. Some special values are 








W (-1,-e7") =-1 
W(-e7') =-l 
W(e) = 
W (0) =0 
(149) W (oo) = 00 
W (—co) =ootin 
w(-§) -4 
W (— In V2)) = — ln(2) 
W (-1,-In(V2)) = -—21n(2) 
We also have the limit 
(150) limg +400 L2 = ami 
and differential 
d W(a,f(2)) ae f(2) 
(151) aV(af2)) = xQaFWaFey) 
as well as the obvious integral 
(152) Jo W (-2@) de = fy -In(a)de = 1 


Let us define, for the sake of brevity, the function 


Win(z) =w (-1,-22) 
(153) In(z) i poo «—3in—In(x)+In(— =) 


Then we have the limits 








lims- Win(x) =0 
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—— R(-In(x)}) — 8(-In(x)) —— |-In(x)| R(W,,(*)) S(r) Mal] 
FIGURE 16. W (-22) = — ln(x) and Wi, (x) = W (-1,-22) 
and 


(155) Im (Win(z)) = | oo 0<r<1 


(156) W(x) =-—In(z) Va ¢ [0,e] 


The root of Re (Wm(x)) is given by 


{x :Re(Win(x)) =O} = 
(157) 


Swine 
Š 
Iw 
S] 


2 
.27441063190284810044... 


~ 
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where the imaginary part of the value at the root of the real part of Wi,(z) is 





(158) =W (-1, 3) 


Oe 
= 14.712388980384689857 ... 


3.1.6. The Lerch Transcendent ®(z,a,v). The Lerch Transcendent [14, 1.11] is defined by 


n 


(159) G2,00) =Po Gin V{|z| < 1}or{|z| = land Re(a) > 1} 


The Riemann zeta function is the special case 


(160) C(s) = ®(1,8,1) = Erco agape 


3.2. Applications of w(x). 


3.2.1. Expansion of y. Consider Euler’s constant y = 0.577215664901533... (33) 
(161) w'(y) = an — bry 


whereupon iteration we see that 


(162) 
n 01 2 3 4 5 6 7 8 9 10 
—an | =| 0 1 48 290 581 1163 2327 13964 7492468716 14984937433 1078915495184 
—bn 1 2 84 504 1008 2016 4032 24192 12980362752 25960725504 1869172236288 
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3.3. Conventions and Symbols. Many of these symbols are from [23, p491]. 





(163) 
i Jl 
R {x : —o0 < x < co} 
R {z : —o0 < x < œ} 
Rt {z:0<4< co} 
R? {z1 ...%q 2-00 < zi < co} 
C {xz + iy: x,y € R} 
Z E S E 
N {0,1,2,3,....} 
N* FDS diac 
H {0,4:neZ} 

f(x) 











f(x) =o(g(a)) lime soo 49 = 0 
f(x) x g(x) fa< 42 <b: {a,b} > 0} 
#A numbers of elements in the finite set A 
|Ala d-dimensional Lebesgue measure(volume) of A C R? 
d(x, A) {min(|x — y|) : y € A} Euclidean distance between x and the nearest point of A 
exp(z) exponential e? = 5? o | 
Res(f(x)) complex residue of f(x) at x = y 
z=y 
|z] floor, the greatest integer < x 
{x} x — |x| , the fractional part of x 
Tz complex conjugate, Re(x) — Im(x) 
Fix} n-th fixed point of the map f(z), n-th solution to f(z) = x 
Pk k-th prime number 
Inp(a) mO) 
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